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Abstract. We study transport of a weakly diffusive pollutant (a passive scalar) 
by thermoconvective flow in a fluid-saturated horizontal porous layer heated from 
below under frozen parametric disorder. In the presence of disorder (random frozen 
inhomogcncities of the heating or of macroscopic properties of the porous matrix), 
spatially localized flow patterns appear below the convective instability threshold of 
the system without disorder. Thermoconvective flows crucially effect the transport 
of a pollutant along the layer, especially when its molecular diffusion is weak. The 
effective (or eddy) diffusivity also allows to observe the transition from a set of localized 
currents to an almost everywhere intense "global" flow. We present results of numerical 
calculation of the effective diffusivity and discuss them in the context of localization 
of fluid currents and the transition to a "global" flow. Our numerical findings are in a 
good agreement with the analytical theory we develop for the limit of a small molecular 
diffusivity and sparse domains of localized currents. Though the results are obtained 
for a specific physical system, they are relevant for a broad variety of fluid dynamical 
systems. 
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The effect of localization in spatially extended linear systems subject to a 
frozen random spacial inhomogeneity of parameters is known as Anderson localization 
(AL). AL has been first discovered and discussed for quantum systems pQ. Later 
on, investigations were extended to diverse branches of classical and semiclassical 
physics: wave optics (e.g., [2J), acoustics (e.g., |3J), etc. The phenomenon has 
been comprehensively studied and well understood mathematically for the Schrodinger 
equation and related mathematical models (e.g., [HO [6]). Also, the role of nonlinearity 
in these models has been addressed in the literature (for instance, destruction of AL by 
nonlinearity [HE]). 

Being well studied for conservative media (or systems) the localization phenomenon 
did not receive a comparable attention for active/dissipative ones as, e.g., in problems 
of thermal convection or reaction-diffusion. The main reason is that the physical 
interpretations of formal solutions to the Schrodinger equation and governing equations 
for active/dissipative media are essentially different and, therefore, the theory of AL 
may be extended to the latter only under certain strong restrictions (this statement is 
discussed in details in the end of the next section). Nevertheless, effects similar to AL can 
be observed in fluid dynamical systems ( [8] ; in [9] the effect of parametric disorder on the 
excitation threshold in one-dimensional Ginzburg-Landau equation has been studied, 
but without attention to localization effects). In this paper, we study an example: the 
problem where localized thermoconvective currents excited under parametric disorder 
crucially influence the process of transport of a passive scalar (e.g., a pollutant). 

The paper is organized as follows. In section [I] we formulate the specific physical 
problem we deal with, introduce the relevant mathematical model, and discuss physical 
background for the problem. Section [2] presents the results of a numerical simulation. 
In section [3] we develop an analytical theory for a certain limit case. Section H] ends the 
paper with conclusions. 

1. Problem formulation and basic equations 

The modified Kuramoto-Sivashinsky equation 

6(x,i) = - (e xxx (x,t)+q(x)6 x (x,t) - (6 x (x,t)f) x (1) 

is relevant for a broad variety of active media where pattern selection occurs. It governs 
two-dimensional (2D) large-scale natural thermal convection in a horizontal fluid layer 
heated from below [101 [11] and holds valid for a turbulent fluid [12] , a binary mixture at 
small Lewis number [13], a porous layer saturated with a fluid [HE], etcU Specifically, 
in the problems mentioned, temperature perturbations 9 are almost uniform along the 
vertical coordinate z and obey ([I]). 

To argue for a general validity of (CQ), let us note the following. Basic laws in 
physics are the conservation ones. This fact quite often results in final equations having 

\ In these fluid dynamical systems, except the turbulent one [12], the plates bounding the layer should 
be nearly thermally insulating for a large-scale convection to arise. 
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the form d t [quantity] + V-[flux of quantity] = 0. With such conservation laws either for 
systems with the sign inversion symmetry of the fields, which is wide spread in physics, 
or for description of a spatiotemporal modulation of an oscillatory mode, the original 
Kuramoto-Sivashinsky equation (e.g., see [15]) should be rewritten in the form (JT]). On 
these grounds, we claim equation (JTJ) to describe pattern formation in a broad variety 
of physical systems. 

In the following we restrict our consideration to the case of convection in a porous 
medium; nevertheless, the most of results may be easily extended to the other physical 
systems mentioned. Equation ([1]) is already dimensionless and below we introduce all 
parameters and variables in appropriate dimensionless forms. 

Recall, the large-scale (or long-wavelength) approximation is identical to the 
approximation of a thin layer and assumes that the characteristic horizontal scales are 
large against the layer height h. For large-scale convection, (21/ 2) h 2 q(x) (cf. (pQ), [H]) 
represents relative deviations of the heating intensity and of the macroscopic properties 
of the porous matrix (porosity, permeability, heat diffusivity, etc.) from the critical 
values for the spatially homogeneous case. Thus, for positive spatially uniform q, 
convection sets up, while for negative q, all the temperature inhomogeneities decay. 
For convection in a porous medium [TJ], the macroscopic fluid velocity field 

v = —e x -—e z , * = —j^-z(h-z)6 x {x,t) = f(z)iJ>{x,t), (2) 

where ip(x,t) = 9 x (x,t) is the stream function amplitude, the reference frame is such 
that z = and z = h are the lower and upper boundaries of the layer, respectively 
(figure [lb). Though the temperature perturbations obey (JTJ for diverse convective 
systems, function f(z), which determines the relation between the flow pattern and the 
temperature perturbation, is specific for each case. 

Though is valid for a large-scale inhomogeneity q(x), which means /il^l/l^l 1 5 
one may set such a hierarchy of small parameters, namely h <C (h\q x \/\q\) 2 -C 1, that a 
frozen random inhomogeneity may be represented by white Gaussian noise £(x): 

q(x) =q + £(x), (£(x)) = 0, (mt(x')) = 2e 2 6(x - x'), 
where e 2 is the disorder intensity and go is the mean supercriticality (i.e. departure from 
the instability threshold of the disorderless system). Numerical simulation reveals only 
steady solutions to establish in ([!]) with such q(x) [8]. 

Let us now discuss some general points related to the physical problem under 
consideration. Obviously, the linearized form of equation (JTJ) in the stationary case, 
i.e., 

-0 xxx (x) - £(x) 9 x (x) = q 6 x (x) , 

is a stationary Schrodinger equation for if) = 6 X with q instead of the state energy 
and —£(2;) instead of the potential. Therefore, like for the Schrodinger equation (e.g., 
see [UEIE]), all the solutions ip(x) to the stationary linearized equation ([I]) are spatially 
localized for arbitrary q ; asymptotically, 

ip(x) oc exp(— j\x\), 
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Figure 1. (a): Establishing steady solutions to |T]) for g indicated in the plot are sets 
of exponentially localized patterns [shown for one and the same realization of random 
inhomogeneity £(x) and e = 1; q(x) is represented by q 7r (x) — tt^ 1 ff^^ q(x')dx'}. 
(b) : The stream lines corresponding to the solutions in graph (a) are plotted for the 
case of convection in a porous layer [cf. |2|]. 



where 7 is the localization exponent. Such a localization can be easily seen for a solution 
to the nonlinear problem ([1]) in figure for q = —2.5. 

One should keep in mind, that, in the quantum Schrodinger equation, localized 
modes are bound states of, e.g., an electron in a disordered media. Even the mutual 
nonlinear interaction of these modes, which appears due to the electron-electron 
interaction and leads to destruction of AL, should be interpreted in the context of 
the specific physical meaning of the quantum wave function. Therefore, the theory 
developed for AL in quantum systems may not be directly extended to active/ dissipative 
media. Indeed, in (CQ), all excited localized modes of the linearized problem do mutually 
interact via nonlinearity in a way where they irreversibly lose their identity (unlike 
solitons in soliton bearing systems, which completely recover after mutual collision). 
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Thus, when the spatial density of excited localized modes is large and these modes form 
an almost everywhere intense flow, localization properties of formal solutions to the 
linearized problem do absolutely not manifest themselves. 

Nevertheless, when excited modes are spatially sparse, solitary exponentially 
localized patterns can be discriminated as reported in [8]. Figure Q] shows sample 
patterns for such a case. One can see that for negative q the spatial density of excited 
modes rapidly decreases as q decreases and the pattern localization becomes more 
pronounced. For a small spatial density of excited modes, one can distinguish all these 
modes and introduce the observable quantifier v of the established steady pattern, which 
measures the spatial density of the domains of excitation of convective flow; fortunately, 
an empiric formula fits perfectly the numerically calculated dependence [8], 



where % = e ' q Q . 

Here we would like to emphasize the fact of existence of convective currents below 
the instability threshold of the disorderless system. These currents may considerably 
and nontrivially affect transport of a pollutant (or other passive scalar), especially when 
its molecular diffusivity is small (for instance, for microorganisms or suspensions the 
diffusion due to Brownian motion is drastically weak against the possible convective 
transport). Transport of a nearly indiffusive passive scalar is the subject of our 
research, as a "substance" which is essentially influenced by these localized currents and, 
thus, provides an opportunity to observe manifestation of disorder-induced phenomena 
discussed in [H]. 

From the viewpoint of mathematical physics, there is one more nontrivial question 
which is worthy to be addressed. In AL an important topological effect takes place; 
while in ID case all the solutions are localized, in higher dimensions spatially unbounded 
solutions appear (e.g., |4j). The modification of ([I]) for the case of inhomogeneity in 
the both horizontal directions, q = q(x,y), (e.g., see [H]) cannot be turned into the 
Schrodinger equation even after linearization in the stationary case. Thus, there are 
no reasons for any topological effects directly analogous to the one mentioned for AL. 
Nevertheless, one may speak of a percolation kind transition, where the domain of an 
intense convective flow becomes globally connected for high enough go- Noteworthy, this 
transition cannot be observed in ID system (CQ) as there is always a finite probability of 
a large domain of negative q(x) where the flow is damped and the domain of an intense 
flow becomes disconnected. Essentially, the flow damped never decays exactly to zero 
and, hence, one needs a formal quantitative criterion for the absence of intense currents 
at a certain point. On the other hand, this transition leads to a crucial enhancement 
of transport of a nearly indiffusive scalar along the layer, and the intensity of this 
transport can be used to detect the transition immediately in the context that arises 
applied interest to it. In this way, one also avoids introducing a formal quantitative 
criterion. Remarkably, in the context of transport of a passive scalar, that is the subject 




(3) 
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of the study we present, the transition from a set of spatially localized currents to an 
almost everywhere intense "global" flow can be observed in ID system ([1]) as well. 

Let us describe the transport of a passive (i.e., not influencing the flow in contrast, 
for instance, to [16]) pollutant by a steady convective flow fl2D- The flux j of the pollutant 
concentration C is 

j = vC - DVC , (4) 

where the first term describes the convective transport and the second one represents the 
molecular diffusion, D is the molecular diffusivity. The establishing steady distributions 
of the pollutant obey 

V-j = 0. (5) 

Equation ((SJ) yields a uniform along z distribution of C (see appendix), 
dC{x) J 



dx „ 21i/j 2 (x) 

D + 



(6) 



2h 2 D 

where J is the constant pollutant flux along the layer. Note, for the other convective 
systems we mentioned above, the result differs only in the factor ahead of ijj 2 /D. 

2. Effective diffusivity 

In this section we introduce and consider the effective diffusivity (for general ideas on 
introducing the effective diffusivity one can consult, e.g., [HI [18]). Let us consider the 
domain x 6 [0, L] with the imposed concentration difference 5C at the ends. Then the 
establishing pollutant flux J is defined by the integral [cf. (jBJ)], 

L 

dx 



5C = -J 



D , 21^ 2 (x) • 



2h 2 D 

For a lengthy domain the specific realization of £(x) becomes insignificant; 
5C=-JL((D + 2 ^Y\=-c--iJL. 



2h 2 D 
Hence, 

5C 

J=-o~, 

i.e. o can be considered as an effective diffusivity. 
The effective diffusivity 

'-(("^)T 

turns into D for vanishing convective flow. For small D the regions of the layer, where 
the flow is damped, ^<1, make large contribution to the mean value appearing in (J7J) 
and diminish a, thus, leading to the locking of the spreading of the pollutant. 



Diffusion of a passive scalar by convective flows under parametric disorder 7 




Figure 2. Dependencies of effective diffusivity a on mean supercriticality go for 
molecular diffusivity D indicated in the plot. The bold black line in the inner plot 
represents the analytical dependence (see section [3]). 



Note, disorder strength e 2 can be excluded from equations by the appropriate 
rescaling of parameters and fields. Thus, the results on the effective diffusivity can 
be comprehensively presented in the terms of D, a, and %: 




Figure [2] provides calculated dependencies of effective diffusivity a on q for moderate 
and small values of molecular diffusivity D. Noteworthy, 

(i) for small D a quite sharp transition of effective diffusivity a between moderate values 
and ones comparable with D occurs near qo = 0, suggesting the transition from an almost 
everywhere intense "global" flow to a set of localized currents to take place; 

(ii) below the instability threshold of the disorderless system, where only sparse localized 
currents are excited, the effective diffusion can be dramatically enhanced by these 
currents; e.g., for D = 10~ 4 , % = —1, the effective diffusivity is increased by one 
order of magnitude compared to the molecular diffusivity. 

Figure [3] shows dependencies of the effective diffusivity on the molecular one for 
nonnegative go- Remarkably, for q > 0.3, the dependencies possess a minimum which 
is in agreement with known general results on interference between turbulent and 
molecular diffusion [19]. The reason is the fact, that for the convective transport the 
molecular diffusion plays a destructive role. Hence, for high qo, where convective flows 



Diffusion of a passive scalar by convective flows under parametric disorder 



8 




10 10 10 

molecular diffusivity D 



10 



Figure 3. Dependencies of the effective diffusivity on the molecular one for 
nonnegative go 



are intense, the enhancement of the convective transport prevails over the weakening of 
the diffusional one as the molecular diffusivity tends to zero; on the contrary, for low qo, 
where convective flows are weak, the decrease of the molecular diffusivity leads to the 
weakening of the transport. 



3. Analytical theory 

Let us now analytically evaluate the effective diffusivity for a small molecular one, 
Del, and sparse domains of excitation of convective flow, v <C 1. We have to 
calculate the average 

P = {{p + ^{x)/D)- 1 ) 

[cf. ([7])] which, due to ergodicity, can be evaluated not only as an average over x for a 
given realization of but also as an average over realizations of at a certain 
point x . Let us set the origin of the x-axis at x . Hence, j3 = ((D + ip 2 (0)/ D)' 1 )^. 

When the two nearest to the origin excitation domains are distant and localized 
near x\ > and x 2 < (see figure @J, 

ip(0) w ^ ie - 7Xl + V> 2 e- 7M , (8) 
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Figure 4. Sketch of two localized flow patterns being nearest to the origin 

where ipxp characterize the amplitude of flows excited around xi j2 . For small D and 
density u, the contribution of the excitation domains to (3 is negligible against the one 
of the regions of a weak flow. Therefore, we do not have to be very accurate with the 
former and may utilize expression ([S]) even for small Xi 



P 



1 



D + ^(0)/D / t 




p{x 1 )p(x 2 ) \ ^ 



D + D- l (4>ie-^ + ^ 2 e-T X2 ) 2 # 

where [^(a^)] is the density of the probability to observe the nearest right [left] 

excitation domain at +x% [— Xq\. For probability distribution P{x\ > x), one finds 
P(xi > x + dx) = P(xi > x) (1 — z/dx), z.e., 4-P(xi > x) = —vP{x\ > x). Hence, 
P(xi > x) = e~ ux , and probability density p(x) = \4--P(xi > x)\ = ue~ ux . With 
regard to averaging over ipx^, it is important that the multiplication of ipi^ by factor 
F is effectively equivalent to the shift of the excitation domain by 7~ 1 lnF which is 
insignificant for F ~ 1 in the limit case we consider. Hence, one can assume ipi^ = ±1 
(the topological difference between ipifa < and ipiip2 > is not to be neglected) and 
rewrite ([9]) as 



P = I / dari / dx 2 z/V^ 1 ^ 







1 



D + D- 1 (e~^ Xl + e"^ 2 ) 2 
1 



+ 



D + D- 1 (e~^ - e-^ 2 ) 2 
These integrals can be evaluated for 1^/7 < 1, and one finds 



■2u_ 

1 



«!>(■!) . (10) 



P \D 

For qo<— 1, one can use the asymptotic expressions for z/ [equation (EJ)] and 

_2 / 1 1 5 , 5 

7 = £ 3 Mg |2 _ - _| 9o | 2 

The latter expression is known from the classical theory of AL (cf. [6]). 

In figure [21 one can see analytic expression ffTO]) to match the numerically evaluated 
a for D = 1CT 4 , go < —1 quite well. With ( flOl) . one can evaluate the convectional 
enhancement of the effective diffusivity below the excitation threshold of the disorderless 
system, and it is given by factor a/D = (2/D) 2u ^ which can be large for small D. 
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4. Conclusion 

Summarizing, we have studied the transport of a pollutant in a horizontal fluid layer 
by spatially localized 2D thermoconvective currents appearing under frozen parametric 
disorder. Though the specific physical system we have considered is a horizontal porous 
layer saturated with a fluid and confined between two nearly thermally insulating plates, 
our results can be trivially extended to a broad variety of fluid dynamical systems (like 
ones studied in pH EED, E21 EG]). We have calculated numerically the dependence of the 
effective diffusivity on the molecular one and the mean supercriticality (see figures [2j [3]). 
In particular, for a nearly indiffusive pollutant (D <C 1), first, we have observed the 
transition from a set of localized flow patterns to an almost everywhere intense "global" 
flow, which results in a soar of the effective diffusivity from values comparable with the 
molecular diffusivity up to moderate ones. Second, we have found convective currents to 
considerably enhance the effective diffusivity even below this transition. For the latter 
effect the analytical theory, which perfectly describes the limit of D <C 1, v <S 1, has 
been developed [equation (flUl) ]. 
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Appendix: Diffusion by stationary flow 

The following derivation of equation ([6]) is performed in the spirit of the standard 
multiscale method (interested readers can consult, e.g., [20j IE]). We consider the 
transport of a pollutant in a layer with boundaries impermeable both for the fluid 
and for the pollutant. In order to derive ([6]), we substitute filtration velocity v from (J2]) 
and pollutant flux j from (j4j) into conservation law (jSJ), and write down 



(the subscripts x and z indicate respective derivatives). The absence of the fluxes of 
the pollutant and the fluid trough the boundary, i.e., j z (z = 0) = j z (z = h) = and 
v z (z = 0) = v z (z = h) = (the superscripts x and z indicate the respective components 
of vectors), results in boundary conditions 



(* Z C) X - (tt x C% - D(C XX + C zz ) = 



(Al) 



z = Q,h: C z = 0. 



(A.2) 



§ Calculation of statistical properties of states of an extensive distributed stochastic system, like the 
one performed in this work, is extremely CPU-time-consuming. 
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We assume D ~ ft, -1 , and use ft as a small parameter of expansion; d z oc ft -1 , 
D = h-'D^, f(z) = h^f^z), C(x,z) = C (x,y) + hCfay) + h 2 C 2 (x,y) + .... 
Then ( 1A.1I) reads 

h-^tf^MCo + hC 1 + ...)] - h- x d z [f^e xx {c + hC 1 + ...)] 

(A.3) 

—ft D_i[Co :XX + Cq >zz + h(Ci tXX + Ci )ZZ ) + ft (C2, xx + C*2 iZZ )+ ...] =0. 
From ( 1A.3I) in the order ft~ 3 : 
—C§,zz = . 

Due to boundary conditions (1A.2I) . Cq )Z = 0, i.e., in the leading order the concentration 
is uniform along z, 

From flA.31) in the order ft~ 2 : 

d x [f-i, z 9 x C ] — d z [f-i9 xx Co] — hD^iCi^ zz = , 

i.e., 

C ltZZ = (ft D^ 1 y 1 f^ hz (z) 6 x (x) r)o, x (x) ■ 
The last equation yields 

C\ = gi(z)r]i(x) + A(x)z + B(x) , 

where 

gx, zz = ^ _1 /-m = 3v / 35(ft - 2z)/h 3 , T] X {x) = Dz\ 9 x {x) r} 0tX (x) . 

Formally, g\ = 3V35(z 2 /2h 2 — z 3 /3h 3 + az + b). Boundary conditions ( 1A.2j) result 
in A(x) = a = 0; B(x) makes a uniform along z contribution to C, like ?7o(x), and, 
therefore, can be treated as a part of r] . Hence, we may claim f C\dz = 0, i.e., 
B(x) = b = 0, and obtain 

( z 2 z 3 \ 
Ci =9i(z) T] x {x) = 3V35 I — - — J ?7i(x) . 

Let us now find the gross pollutant flux ft J through a vertical cross-section of the 
layer. The integral of j x [equation (jlj)] over z is 

h h h 

h,J= J fdz = J (v x C - DC x )dz = J (V Z C - DC x )dz 



h h 

= J [f z o x c - Dc x )dz = J (f z e xVo + h f z9l e xVl - d Vo>x + o(h 3 ))dz 



h 

~ h9 x r]i j fz9idz ~hDr] 0:X = -ft (-^^(9 x ) 2 r] , x + Dr] ^ . 



o 



As the pollutant is not accumulated anywhere, flux J should be constant along the layer. 
Hence, we find equation ([6]) providing the relation between the concentration field and 
the flux. 
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